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We present results of an extension of the market model introduced by 
Bornholdt to high dimensions. Three and four dimensions are shown to behave 
similar to two, for suitable parameters. 
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1. Introduction 

To construct a reasonable dynamical financial market model consistent with the 
observed properties in real market is a central theme. Many models have been 
proposed along this line. Specifically application of concepts and techniques of sta- 
tistical mechanics are giving deeper insight into the understanding of the complex 
behavior of a market. Fluctuations of prices of commodities, stocks and foreign 
exchange rates are examples of these. Reproduction of the distribution of returns 
based on the empirical data are a strong criterion to select models. It is important 
to set a communication structure for the model of financial market. Various models 
are proposed to iniitate the real market commrmjijiiljiin Among them we mention 
the minority gameld, percolationnB , Ising modelfflu'ErEltl and so forth. In this paper 
we focus on a spin model that seems to be appropriate to capture the movements 
of constituent agents in a simplified manner. In previous work, lorU has modeled 
three possible states of each agent's decision as spin-1 model and its decision mak- 
ing is affected by two different kind of noises and their nearest neighbors.cl However 
in this model the influence of market prices on each trader decision is not incor- 
porated. Bornholdt recently formulated these influences with two different scales 
of local interaction and global one in two dimension.! A metastable phase appears 



*On leave from Department of Applied Physics, Faculty of Science, Tokyo Institute of Technology, 
Oh-okayama, Meguro-ku Tokyo, 152-8551, Japan 
E-mail: tyamano&mikan. ap.titech. ac.jp 



T. Yamano 



at randomly frozen finite magnetization which correspond to a bubble-Hke state in 
terms of economics.S 

We assiune that the traders hve on a d-dimcnsional lattice. This totally contrasts 
the casaj where the super-spins do not need any geometry. In this paper we present 
the results of higher dimensions in the Bornholdt model. This can be motivated 
from the fact that one of the most important characters of the market evolution is 
the heterogeneous structure of agents. The heterogeneity arises from not only the 
geometrical distance of agents but also the amount of knowledge that each agents 
has, accessibility of market information prevailing, preferences, and processing skills 
for investments etc. 



2. The Bornholdt Model 

In the model, the spin variables are either +1 or —1, which allows each agent 
to decide actions of buying or selling at each time step t. Thus the magne- 
tizetion M{t) = J2j=iSj{t)/N can be interpreted as a measure of priceE We 
identify the logarithm of the absolute value of the magnetization as returns, i.e. 
ret{t) = ln|M(i)| — In |M(i — 1)|. The use of linear returns (or simply relative 
price chanffc) is also possible and in fact the main properties are also seen for this 
definition. tj Each spin is updated by the following heat bath dynamics, 

^ r +1 withp= [l + exp(-2/3/i,(t))]-i 

1 ~1 with 1 — p ' 

where /3 is the inverse temperature and hi{t) is a time-dependent local fieldjthat each 
spin feels. This corresponds to a signaLthat each agent i receives at time m and also 
to "individual bias" of the i-th agent .Q The temperature is introduced in a totally 
parallel way to the Ising spin model in magnetism. This fictitious temperature has 
been used in the previous works. BlI 

Bornholdt introduced the following form of the local field, 

N N 

The above form can be considered incorporating two contradictory movements of 
constituent agents in a real market: herd behavior and preference of minority. J.^ 
denotes the strength of the communication between the agents i and j. Further- 
more Jij is set to a constant J for nearest neighbor agents and zero for all the other 
pairs. This term corresponds to the "disagreement function" of the i-th agentB a 
is a strength of the coupling to the magnetization and assumed to be positive. In a 
spin language, the first term prefers the ferromagnetic state. The second term, on 
the other hand, tends to encourage a spin flip when magnetizations becomes large.! 
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The more complicated form of the local field is given asE 

N ^ N 



N N 

h^ (t) = J,, S, it) - aa it) - J2 Si it) , (3) 

j=i 

meaning we have two spins on each site at each time: demand spin Si{t) and strategy 
spin Ci{t) of an agent i. Consideration in a market context gives an interpretation 
of strategy spin asHI 

„ , , J +1 anti-ferro (fundamentalist) 

" \ -1 ferro (chartist) ^ > 

This reminds us of the Lux-Marchesi modelO, where agents are divided into these 
two groups and evolve with interacting mutually and moreover with changing their 
strategies. We flip the strategy spin at next time step: Ci{t + 1) = —Ci{t) if the 
Si{t) is antiparallel to the global coupling to the magnetization i.e., 



1 ^ 

SMC^{t)-Y,S,{t)<G (5) 



3 = 



is satisfied at each time because a positive energy contribution (risky or preference 
being minority) to the dynamics encourages these fundamentalists to switch their 
strategies to chartists and vice versa. These updates are done after updating Si{t). 

3. Simulations and Results 

We performed the simulation with the reduced version of the local field eq.(^ 
except for dynamical evolutions of the fundamentalists-chartists behavior. The 
same behavior is observed either with eq. (||) or with eq. (|^) . The size of the hyper 
cubic lattice was set at L = 101, 21, and 7 for the dimension d = 2, 3 and 4 
respectively. We have chosen the temperature (in units of J/ks) T = 1.5, 4.0 and 
6.0 for the dimension c? = 2, 3 and 4 respectively. Note that these temperature is 
below the critical temperature Tc = 2.269(2(i), 4.511(3(i) and 6.680(4d) in the case of 
a = . This selection of parameters allows us to observe the intermittent behavior 
in temporal evolution of ret{t) or clustered volatilities. Fig.l shows the distribution 
of returns for three different dimensions, whose shapes are non-Gaussian. The 



volatility of the returns or the genep 
measured with the following quantityt 



ized cumulative absolute returnS can be 




V{t) = - J2 H{t')\\ (6) 
n ^ — ' ' ' 



t'=t 



where n is the number of data sampled and 7 is a real exponent. We set n = 1 
and 7 = 1, which leads to the cumulative distribution of absolute returns (Fig. 2). 
The measured exponents of these curves are —1.15, —1.41 and —1.50 respectively 
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using the intermediate region of 8 . . . 150, 15 . . . 150 and 40 . . . 150 for fitting, where 
power-law scaling can be seen. A volatility clustering can be seen as a slow decay 
of the autocorrelation function of absolute returns in Fig. 3 on a log-log plot. Fig4. 
is the replot of the same data obtained in our simulation. The slope of the each 
curves is determined as — 1.85 x lO^'^, —1.80 x 10~^ and —7.90 x 10"^ with the fitting 
time 350 . . . 1000(2d), 40 . . . 140(3d) and 5 . . . 35(4d) respectively, where the accurate 
value changes as the fitting regions vary. Fig. 4 shows the ratio of the difference 
between the number of fundamentalists and that of chartists to the total number 
of agents with eq(^ in the 4d case. The decrease in the ratio of fundamentalists 
corresponds to high variance in changes of returns. 



4. Discussion and Summary 

In this paper, we have performed simulations of the recently proposed model by 
Bornholdt in hyper cubic lattices. We confirmed that there is no intermittent be- 
havior in return(magnetization) when we employ only the second term of eq.(^). 
This means that market intermittency emerge as a consequence of a frustration 
across two different scales: local connection with neighbours and global coupling 
with marketB The effects of dimensionality was confirmed, which exhibits similar 
behavior for three different cases. Moreover we confirmed that the similar inter- 
mittent behavior in returns appear when we drop out Si in the second term of the 
local field (eq.(Q)). In this case, however, the fictitious temperature of the mar- 
ket becomes higher than that of the original one(r = 7.5, 9.3 and 11.0 in the 2d, 
Sd, and 4(i case respectively with the same other parameter values as in the text). 
However their temperatures exceed the critical ones. The exponent of cumulatiii 
distribution of returns is one of the ongoing discussion points for many models .ut 
The recent measurement of the exponent based on a empirical data (S&P 500)t 
says the asymptotic slope of close to —3. A percolation model provides —2.9 with 
201 X 201 latticeO and —4 in high dimensions. In this sense, percolation model gives 
reasonable distribution. On the one hand, the modelof Biham et al. provided the 
range of —2.5 ... — 3.5 for the tails of the distributionjlj which depends on the num- 
ber of agents and a factor characterizing the market dynamics. The slope (moreover 
the shape) of the cumulative distribution is strongly affected by the temperature in 
the present model. We confirmed the effect in the 3d case by fixing the other param- 
eters. This means the power-law scaling in intermediate regions is very sensitive to 
the temperature. More specifically, cumulative distributions have a shoulder when 
the temperature is not the tuned value. The intermittency becomes more sparse as 
the temperature descends. The same can be seen as the lattice size becomes large. 
As for autocorrelation function, the empirical data are_well fitted with the exponent 
—0.3 according to Gppikrishnan et al.EJ and Liu et al.EJ. The exponent determined 
from minority gamea is —0.64. In the present model, the tails of the autocorrelation 
function do not behave like power laws but as exponential up to four dimensions. 
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Fig. 1. Distribution of the return defined as the difference of the logarithm of the absolute mag- 
netization. Parameters J = 1.0 and a = 4.0 are common for all dimensions. The symbols 0, + 
and X correspond to {d,T) = (2, 1.5), (3,4.0), and (4,6.0) respectively. 
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Fig. 2. Cumulative distribution of absolute returns multiplied by 100. For all lines, ./ = 1.0, a = 4.0. 
The solid, broken and dotted line correspond to 2d with T = 1.5, 3d with T = 3.8 and 4d with 
T = 6.0 respectively. For better statistics, we have summed over lO'^ Monte Carlo steps. 
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Fig. 4. Autocorrelation function of absolute returns in semi-log scale replotted. The fitted straight 
line to the 2d data has exponent —1.85 x 10~^. 
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Fig. 5. The change of the market constituents in 4d case (J = 1.0, a = 4.0 and T = 6.0) at a 
period of 4000. 
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